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Starting with empirical tight-binding band structures, the branch-point (BP) energies and resulting valence 
band offsets (VBOs) for the zincblende phase of InN, GaN and AIN are calculated from their k-averaged 
midgap energy. Furthermore, the directional dependence of the BPs of GaN and AIN is discussed using 
the Green's function method of Tersoff. We then show how to obtain the BPs for binary semiconductor 
alloys within a band-diagonal representation of the coherent potential approximation (CPA) and apply this 
method to cubic AIGaN alloys. The resulting band offsets show good agreement to available experimental 
and theoretical data from the literature. Our results can be used to determine the band alignment in isovalent 
heterostructures involving pure cubic Ill-nitrides or AIGaN alloys for arbitrary concentrations. 
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I. INTRODUCTION 

Due to their wide range of potential application in op- 
toelectronic and electronic devices, the III- nitride based 
semiconductor systems InN, GaN and AIN and their 
ternary and quarternary alloys have attracted a great 
deal of interest for many years now. ' They can crystal- 
lize in the thermodynamically stable hexagonal configu- 
ration with a wurtzite crystal structure and in a cubic 
modification with a zincblende structure. Due to the low 
symmetry of the wurtzite system, hexagonal nitrides do 
not only show strain-induced, but also spontaneous po- 
larization. The latter effect is absent in the zincblende 
systems: it is for example possible to eliminate the built- 
in fields in quantum dots by the growth of cubic III- 
nitride along the nonpolar (001) direction, e. g. by using 
GaN as dot and AIN as barrier material. " Furthermore, 
quantum well and superlattice structures based on cu- 
bic nitrides can be grown with good structural quality. 
AlGaN/GaN heterostructures can e.g. be used to cre- 
ate a two-dimensional electron gas without polarization 
effects. 

For the design as well as the theoretical simulation of 
such semiconductor devices, the knowledge of the respec- 
tive valence band offsets (VBOs) across the junction of 
two semiconductor materials A and B is of crucial in- 
terest. In low-dimensional structures like superlattices, 
quantum wells or quantum dots, the confinement poten- 
tial is realized by a A/B hetero junction that results in a 
suitable conduction band (CB) and valence band (VB) 
offset. Depending on the values of the respective band 
gaps and the resulting VBO, the electrons and holes are 
either located in the same ("type I" alignment) or in dif- 
ferent spatial regions ("type II" alignment). 

Furthermore, the VBO does not only determine the 
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geometry of the confinement potential and thus the elec- 
tronic properties of low-dimensional systems. It is also 
relevant for the simulation of alloyed systems of the type 
Aa;Bi_a;, because the constituents have to be aligned 
along a common energy reference scale to begin with. 
The concentration dependence of e. g. the band gap will 
therefore directly depend on the energetic offset. In sys- 
tems like the above mentioned AlGaN/GaN, we even find 
a combination of these two effects. 

Although many different experimental and theoretical 
methods are employed to determine VBOs, the disparity 
in the literature values is often large. On the theoreti- 
cal side, computationally costful ab initio methods^' ' are 
ever-improving, but often still lack quantitative agree- 
ment with the experiment and do not really provide the- 
oretical insight into the underlying formation mechanism. 
A different approach is presented by models that predict 
band discontinuities by alignment to a common energy 
reference like vacuum levels, or different kinds of 
charge neutrality levels. '^"'"^^ This determination of tran- 
sitive "natural" VBOs represents a rather macroscopical 
approach that can give insight in the formation mech- 
anism, but neglects small-scale effects like lattice mis- 
match or band bending on small length scales at the 
boundary. 

In the present work, we use two approaches from the 
literature to determine the charge neutrality level, also 
known as branch point (BP), of the continuum of in- 
terface states of cubic InN, GaN and AIN with band 
structures from a flexible empirical tight-binding model 
(ETBM). We then calculate the resulting VBOs between 
the Ill-nitrides and discuss their directional dependence 
with the help of a Green's function method from the 
literature. The experimentally observed directional in- 
dependence of the cubic GaN/AlN VBO can nicely be 
explained in terms of our results for the BPs. 

As our results for the pure systems agree well with 
literature data (where available), we use the coher- 
ent potential approximation (CPA) to calculate the 
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concentration-dependent BP position for disordered cu- 
bic AlGaN alloys and compare it to results from the 
simpler, but often applied virtual crystal approximation 
(VGA) . While the VGA fails to reproduce the concentra- 
tion dependence of the BP, the GPA results yield good 
agreement with experimental data for Alo.3Gao.7N/GaN 
interfaces. Our calculated GPA BPs can be used to esti- 
mate the band alignment in interfaces with AlGaN alloys 
for arbitrary concentrations. 



II. THEORY 

A. Branch-point energy and band alignment 

The simplest linear models for the determination of 
VBOs across common interfaces between semiconductors 
use the alignment along a common reference level. An in- 
tuitive approach is the use of the vacuum level, i. e. the 
comparison of the electron affinities of the constituent 
materials.^" However, the presence of interfaces and sur- 
faces will influence the electronic properties, so that the 
use of values for infinite bulk crystals does not give satis- 
factory results for semiconductors. A large body of work 
deals with the refinement of this approach by establish- 
ing suitable vacuum reference levels for model solids (see 
the excellent review part of Ref. for further reference). 

Another approach that renders good results is the use 
of a common charge neutrality level. It is for exam- 
ple possible to align a broad range of materials along 
neutrality levels of foreign atoms like interstitial hydro- 
gen.^'' In a similar manner, it is well established to em- 
ploy the charge neutrality level of the continuum of inter- 
face states. These localized states develop from itinerant 
Bloch states, but are characterized by a complex wave 
vector k, where Im(k) is normal to the plane that breaks 
the translational symmetry. Gonsequently, they decay 
exponentially into the other material and carry charge 
density across the interface, which induces an interface 
dipole."^" The charge neutrality level of the spectrum of 
interface states is commonly called the branch point en- 
ergy ii^BP and is basically the crossover level between a 
donorlike and acceptor like character. The Fermi level 
will be pinned near the BP and the net charge density 
transfer due to interface states will depend on the sign 
and magnitude of Ep — E^p. The (hypothetical) exact 
alignment of the BPs of two materials with a common 
interface corresponds to a band lineup where all dipole 
contributions cancel each other out. 

In Ref. Tersoff noted that Ebp can be calculated 
from the bulk properties of a semiconductor by the use 
of the cell-averaged real-space Green's function for the 
propagation along a lattice vector R. This function is 
given as 
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with En(k.) as the dispersion of the n-th band. We de- 
note this Green's function with the subscript "int" to 
avoid confusion with the later introduced GPA Green's 
function, although the definition (1) is per se independent 
of any interface. 

Because donor states predominantly take their spectral 
weight from the VB, while acceptor states likewise orig- 
inate from the GB, E^p can be identified as the energy 
where the GBs and VBs equally contribute to G^^{E) 
when the interface normal is chosen parallel to the direct 
lattice vector R. As Tersoff pointed out, it is sufficient to 
integrate over the three-dimensional bulk Brillouin zone 
(BZ) instead of separating ky and kj^, as the choice of a 
large enough R automatically projects out the relevant 
contribution to G^^{E). This method is not applicable 
when the BP lies outside the band gap, however. 

When the orientational dependence is small, the above 
mentioned approach can substantially be simplified by 
the use of interface-averaged approximations. In an at- 
tempt to generalize previous approaches,''''-"^'-'"' Schleife 
et al. calculated the BP as a BZ average of the midgap 
energy. Reformulated in terms of the band-resolved 
densities of states (DOS) (ffve) ^ number of 
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This formula also holds when the BP lies in a band, as 
experimentally known for some small-gap materials like 
InAs or InN. As no interface direction breaks the 
symmetry in Eq. (2), the BP can in principle also be ob- 
tained from the energies at the Baldcreschi point,''"' see 
e. g. Ref. ■->2. 

It should be noted that both methods do not only re- 
quire an appropriate input band structure; in practice, 
the choice of a proper subset of conduction and valence 
bands around the band gap region is also necessary. 

Strictly speaking, the knowledge of the energetic posi- 
tion of the BP is not sufficient to determine the VBO. In 
real hetero junctions between two materials A and B, the 
local charge neutrality condition will always be violated 
to some extent. The difference in the electron affinities 
results in a net dipole, which itself will be screened by 
a characteristic dielectric response e. If the BP energy 
i?BP of each material is measured with respect to its VB 
edge, the VBO AE^ is obtained as 
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The magnitude of the interface dipole contribution -Edip 
can be estimated following the work of Monch.- ' With 
Xa and Xb as electronegativity values from the Miedema 
scale, we have 



Edip = Dx(Xb — Xa), 



(4) 
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with the phenomenological slope parameter 



Here, = 0.86 eV/Miedema-unit connects the work 
functions and the electronegativity scale, while £00 is 
the high-frequency dielectric constant of the semiconduc- 
tor. The Miedema electronegativity values were origi- 
nally defined for elementary solids. Suitable values for 
binary compounds semiconductors can approximately be 
obtained from the geometric mean of the constituents' 
values, e. g. using the data from Table A. 4 of Ref. 
Note that Eq. (4) has originally been deduced for metal- 
semiconductor contacts ' but was then found to be ap- 
plicable to semiconductor interfaces, too. 

B. Empirical tight-binding model 

The importance of the choice of the band structure 
method in BP calculations has already been discussed 
in Ref. . It has been pointed out that a suitably 
parametrized empirical tight-binding model will allow for 
a flexible flt to a large parameter space; the input param- 
eters can then be either chosen from experiments or from 
results of more sophisticated calculation schemes. They 
can as well deliver a more or less realistic dispersion over 
the whole BZ, including the Van Hove singularities at 
the zone boundaries, which give relevant contributions. 
Furthermore, due to the low computational cost, dense 
wave vector samples samples on large k-space regions are 
effortlessly available, which is crucial to the convergence 
behaviour of the Green function, Eq. (1). 

Throughout this paper we will use a tight-binding 
model based on the work of Loehr.^' It makes use of 
a Wannier-like sp^ basis per spin direction and lattice 
site and can thus replicate eight doubly degenerate bands 
for crystals with zincblende structure. The tight-binding 
matrix elements of the bulk Hamiltonian iJ^)""^ are 

iJ^i^' = (Ra|i/'^""^|RV), (6) 

where R runs over the sites of the underlying fee Bravais 
lattice and a is the orbital index. The band structure 
i?„(k) is obtained by the diagonalization of the matrix 

a' R 

for each k e 1. BZ. 

By restriction of the matrix elements to the second- 
nearest neighbor shell, it is possible to fit the band struc- 
ture to 11 material-specific parameters, which are listed 
in Tab. I (the reader is referred to Refs. 26, 37, and 38 
for a more detailled discussion). The fact that all bands 
are also fitted to the X-point energies makes it especially 
suitable for indirect materials, of which the properties are 
often not yet unambiguously established. One notable 



example is AIN, which has a direct band gap in its hexag- 
onal modification, while an indirect T-X band gap is 
commonly " but not always ' assumed for the zincblende 
phase. 

Note that a modification of the model such as in Ref. 26 
is not necessary here, as the spin-orbit interaction is com- 
paratively weak in the Ill-nitrides, which leads to spin 
splittings two orders of magnitude smaller than the band 
gap. 

In order to ensure the necessary consistency among 
the input values (as e. g. the energetic distance of bands 
influences their curvature due to band-mixing effects) 
and simultaneously avoid the common band gap prob- 
lem in ab initio calculations, we use a comprehensive set 
of material parameters from Fritsch et al. from empirical 
pseudopotential calculations (Ref. for GaN and AIN, 
Ref. 40 for InN). The input parameters for the pseudopo- 
tentials stem from various sources which are listed in the 
mentioned references. In case of InN, we additionally use 
recent OEPx(cLDA) + GoWq band curvature parame- 
ters TOc,7i,72,73 from Rinke et al., ' as their ab initio 
band gap for zb-InN does only deviate slightly from ex- 
perimental results. Furthermore, we followed the advice 
of Vurgaftman and Meyer and corrected the AIN gap 
slightly upwards, '" as the original value from Thompson 
et al. has originally been determined for room temper- 
ature, ' ' while the remainder of the input parameters is 
assumed to be valid for low temperatures. 

The resulting ETBM band structures and band- 
resolved densities of states for zb-GaN and zb-AlN are 
depicted in Fig. 1, together with the first moments /cen- 
ters of gravity of the bands (the InN results are not added 
for the sake of clarity). 

C. Coherent potential approximation in band-diagonal 
representation 

The CFA"''"^' is a well-established self-consistent 
Green's function method for the treatment of substitu- 
tionally disordered systems. The CPA itself and its broad 
scope of application is excellently described in detail in 
Ref. 4N. For the sake of completeness, we will here only 
repeat the most important features. 

In the standard form of the CPA, the TB Hamiltonian 
H of the Ax^i-x alloy is seperated into a site-diagonal 
part V and an off-diagonal part W, such that H = V+W. 
Here, only V is considered to be site-dependent, depend- 
ing on the species on the site R. The infiuence of the 
site-diagonal disorder is absorbed into the self-energy op- 
erator 

j:{z) = H,siz)-W, (8) 

where z is the complex energy and Hcs{z) is an effective 
Hamiltonian; it is defined such that the configurational 
average (. . .) over the resolvent of H equals the resolvent 
of Hctfiz): 

{{zl - H)-') = [zl - H,ff{z)r\ (9) 
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Table I. Input material parameters, taken from Refs. 40, 42, and 4-5. See text for further details. 



Parameter 


Description 




zb-AlN 


zb-GaN 


zb-InN 


a 


lattice constant 


(A) 


4.38 


4.52 


4.98 


Aso 


spin-orbit splitting 


(eV) 


0.019 


0.017 


0.006 


71 


Luttinger parameter 




1.85 


2.89 


6.82 


72 


Luttinger parameter 




0.43 


0.85 


2.81 


73 


Luttinger parameter 




0.74 


1.20 


3.12 




CB effective mass 


1^0 j 


\J.16 


U.i4 


U.Uo 


T\ {Tl) 


CB energy 


(eV) 


5.84 


3.31 


0.59 


ris 


HH/LH VB energy 


(eV) 













CB energy 


(eV) 


5.44 


4.43 


4.76 


XI 


HH/LH VB energy 


(eV) 


-2.32 


-2.46 


-1.48 


XI 


split-off VB energy 


(eV) 


-5.39 


-6.29 


-4.80 




When a representation is chosen, this leads to a set of self- 
consistency conditions for the self energies, the so-called 
CPA equations. In the CPA, the self-energy operator is 
diagonal in every representation. Additionally, its matrix 
elements are neither dependent on k nor on R. 

The combination of the CPA with the here employed 
multiband ETBM has already been presented in Ref. 
together with a critical discussion of the extent of its va- 
lidity. This reference also describes in detail how the CPA 
self-energy and Green's function can properly be trans- 
formed from the localized sp^ basis to a band-diagonal 
representation. The latter, in which H^s is diagonal in 
k, is crucial for the calculation of BPs for the alloyed 
system, as one has to restrict the calculation to a subset 
of bands and the most natural choice of band index n is 
their energetic hierarchy. 

In the here employed form, the CPA only maps the 
substitutional disorder in the first moments / centers of 
gravity i^^yg of the constituents on the self energy 



^Ti,A/B 



(z), (10) 



where (?^/b again the band-resolved DOS (normed to 
unity). The off-diagonal disorder effects are simulated 
in the VCA, i. e. interpolated linearly with the concen- 
tration X. The self-consistent iteration scheme for the 
CPA equations as described in detail in Sec. II of Ref. !8 
finally leaves us with the CPA Green's function. Using 
its k-representation Cqpa, we can then extract the CPA 
one-particle spectral function Sgp^ and DOS g^pj^{E): 

S^^AiE) = -- hm ImGS^A(^ + ^S), (11) 

TT i5\,0 

9CPA (^) = - 4r 1^0 E ^'cPA + • (12) 

Here, (5 € M and iVk is the number of wave vector values 
in the irreducible wedge. 
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Although the CPA can simulate many effects that go 
beyond the scope of static mean-field treatments like the 
VGA (such as spectral broadening due to finite lifetimes 
or band gap bowing), it also has several limitations. Due 
to the translational invariance of i?off, the CPA cannot 
simulate local strain patterns. It will therefore only yield 
reliable results if the difference of the lattice constants 
is rather small. Besides this constraint, the CPA can be 
shown ''^ to yield very good results in two limit cases: 

(1) The weak scattering limit, where the difference of 
the first moments of the A and B bands is smaller than 
the respective bandwidths. 

(2) The split band limit, where these first moments are 
either far apart and/or the bandwidths are small enough 
for the bands of A and B to not overlap. 

D. Calculation of branch points In the coherent potential 
approximation 

The approximated branch point energies for alloys can 
also be calculated from the ETBM+CPA results. We 
first define the complex quasi-particle band structure 

i?CPA(k) = (kn|77efi'|kn) 

= (kn|Ty|kn)+I]"(z==£;,'rPA), (13) 

where |kn) are Bloch states which diagonalize Hcs- Now 
the interface Green's function for the alloy can obviously 
be defined by using i?^^^(k) as dispersion in Eq. (1). As 
the CPA Green's function itself is given by 

GSpa(^) = (kn| [zt - H,s{z)r' |kn) , (14) 

the interface Green's function can directly be obtained 
by means of a discrete Fourier transformation: 

GS,(i?)=^e^'^-^[i?-£;r^k)]-^ 

nk 

= J2e^-^G^^l^iE). (15) 

nk 

Like in the pure case, the BP has to lie in a real gap of 
the quasi-particle band structure. As Ggp^^ has no poles 
in this region, the calculation can be restricted to the real 
part of the energy axis. 

When the BZ average approach from Eq. (2) is used, 
the band-rcsolved densities of states can analogously be 
replaced by their CPA counterparts from Eq. (12). 

III. RESULTS 

A. Pure zincbiende Ill-nitrides 

In order to calculate the BP energies from the TB band 
structures, we follow the recommendations of Schleife et 
al. ^ ' and use one CB and the two uppermost VBs per 
spin direction. As Fig. 1 suggests, this band subset can 
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Figure 2. (Color online) Real part of GjJt ior AIN for the 
(001), (110) and (111) direction. 



be justified by the much larger dispersion of the low- 
est VB. It should however be kept in mind that this set 
is obviously not suitable for several II- VI materials, ""^ 
and the choice of bands certainly remains a major source 
of uncertainties in the BP calculation. For that reason, 
we carefully recommend the same uncertainty range of 
0.2 eV as given by Schleife et al. 

For the BZ average method, convergency of the results 
for Eqp is easily achieved by using a moderately dense 
k-grid on the irreducible wedge, so that a reduction to 
the Baldereschi point does not render any practical ad- 
vantage. Convergency within 10 meV could be obtained 
with a equidistant resolution of about 100-200 k-values 
between F and X. 

The interface Green's function however is known to be 
extremely difficult to converge.^''-'' In practice, it turned 
out that very dense wave vector samples on the full BZ 
are necessary. Figure 2 exemplarily shows the real part of 
^St fo'^ ^O'" three different interface orientations R oc 
(001), (110) and (111). As explained in Sec. II A, the BP 
can be obtained by the zero-crossing of this curve, which 
should arise for large enough multiples of R. For GaN 
and AIN, the intersections could be brought to conver- 
gency within 20 meV for the (001) and (110) direction 
when ~250-300 k-values were used between F and X. 
However, satisfactory convergency for the (111) direction 
could not be achieved for either material. 

In case of InN, the results from the BZ average method 
render a BP energy in the CB region, so that the Green's 
function method is not applicable. This is also in agree- 
ment with theoretical and experimental results from the 
literature.-'^ 

The results for the BP energies Ebp for the cubic III- 
nitrides InN, GaN and AIN can be found in Tab. II, to- 
gether with several theoretical and experimental results 
from the literature. Further information on the respec- 
tive source is provided in the footnotes of the table. 
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Table II. Branch point energies Ebp for cubic Ill-nitrides in eV. 



Material 








-Ebp 








R oc (001) 


R oc (110) 


BZ average 


Literature 


zb-AIN 


5.44 


3.38 


3.64 


3.28 


2.97f 3.19±0.18'' 


zb-GaN 


3.31 


2.62 


2.86 


2.50 


2.37? 2.34±0.09;' 2.31"^ 


zb-InN 


0.59 






1.65 


1.51f 1.59±0.26^ 1.50f 1.38±0.1'= 



" ETBM + BZ average calculation, Ref. no explicit specification on bulk parameters. 

^ Empirical value from fit to experimental band offsets, Ref. . 

" Quasiparticle BZ average calculation with Eg = 3.28 eV (IfSEOS+Go Wo), Ref. 17. 

Quasiparticle BZ average calculation with Eg = 0.47 eV (HSE03+Go Wo), Ref. 17. 
" XPS measurement on (001) surface at room temperature, Ref. 49. 



However, it should be noted that the energy gap of the 
theoretical literature values does in general not coincide 
with our input parameter set. Furthermore, no direction- 
dependent BP energies are available to our knowledge. 
Still, our calculated values agree reasonably well with the 
external data. For the sake of completeness, it should be 
noted that Ref. " contains further theoretical values for 
the cubic Ill-nitrides that were obtained from slightly dif- 
ferent ab initio band energies, but also show good agree- 
ment. 

Our results also show that the position of Ebp carries 
a clear-cut directional dependency, where 

£;bp(BZ average) < -BbpIOOI) < £;bp(110) 

for both GaN and AIN. Additionally, the trend of the 
available data for the (111) direction (see again Fig. 2) 
also suggests that 

Sbp(III) > Sbp(IIO) > £;bp(001) 

might hold. Whether the "real" directional average of the 
BP in these materials lies energetically below the (001) 
and (110) values or the BZ average method just gives 
a slightly too low E-qp, cannot ultimately be concluded 
from our calculations, as reliable results for further di- 
rections would be necessary. 

With knowledge of the BP energies, the band align- 
ment can now be calculated from Eqs. (3)- (5). For the 
contribution of the interface dipole, the electronegativ- 
ity values from Ref. 34 have been used. As the localized 
interface states exponentially decay into the energy bar- 
rier caused by the VBO, we used the dielectric constant 
eoo(InN) = 7.0 '- for both the InN/AlN and InN/GaN 
junction and eoo(GaN) = 5.3 for GaN/AlN. The re- 
sults for the BP difference, the dipole contribution and 
the resulting VBO can be found in Tab. Ill, together with 
theoretical and experimental values from the literature. 

For all interfaces, the contribution of the dipole correc- 
tion is small, so that an approximate alignment off these 
three systems along the BP energy is well justified. This 
is in agreement with general considerations in the liter- 
ature on junctions between cubic binary semiconductor 
compounds.""'-" 

For the GaN/AlN interface, the VBO is practically in- 
dependent of the interface orientation. This fact might 



seem strange at first glance, as the zincblende (001) inter- 
face is polar, while the (110) is not. However, this result 
is in compliance with various experimental and theoreti- 
cal findings (see e. g. Ref. 39) and has also been deduced 
by Lambrecht and Segall~ ' from an analysis of the bond- 
ing geometry and polarity at isovalent zincblende inter- 
faces. Our calculations show that this independence can 
also be deduced from band-structure related properties, 
as the difference in the (001) and (110) BP energies is 
almost identical (see again Tab. II). 

Notably, the BZ average method also gives the same 
BP difference and thus VBO, which underlines its suit- 
ability for material systems with weak directional depen- 
dence. Although no directionally dependent calculation 
could be performed for InN, the detailed analysis of Lam- 
brecht and Segall indicates that the same behavior will 
also hold for InN. 

From a quantitative point of view, the VBOs agree 
quite well with available data from the literature. Our 
value of A£;vB (GaN/AlN) sa 0.8 eV agrees with the rec- 
ommendation from the topical review of Vurgaftman and 
Meyer and also with the comprehensive experimental 
analysis of Monch. The theoretical value of 0.6 eV from 
Monch also stems from ETBM calculations, but with the 
use of atomic term values and hopping matrix elements 
from 1972 and 1981, respectively.""'''" The smallest lit- 
erature value of 0.5 ± 0.1 eV has been determined from 
lattice-matched quantum wells on 3C-SiC substrate at 
room temperature. For both InN/AlN and InN/GaN, 
our results also agree well with the literature, only the 
review by Vurgaftman et al. recommends smaller values. 
This could be explained by the fact that their publica- 
tion is from 2003, when high-quality wurtzite as well as 
zincblende InN samples did not yet exist for a long time. 

The overall resulting band lineup as recommended by 
us for the cubic phases of InN, GaN and AIN is depicted 
in Fig. 3. 

B. AIGaN alloys with zincblende phase 

In contrast to InGaN and InAlN, the lattice mismatch 
between the constituents in cubic AIGaN is relatively 
small with 3%, which suggests that the influence of local 
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Table III. Resulting valence band offsets AEvsiA/ B) for cubic nitrides in eV. 



Interface A/B Eip - Eip 




Dipolc 




AEysiA/B) 


Roc (001) Roc (110) 


BZ average 




This work 


Literature 


GaN/AlN 0.76 0.78 


0.78 


+ 0.02 


0.78, 0.80 


0.5±0.1f 0.6^ 0.8'; 0.85"* 


InN/AlN — — 


1.63 


+ 0.03 


1.66 


1.15^ 1.46^ 1.60=^ 


InN/GaN — — 


0.85 


+ 0.02 


0.87 


0.3" 0.75f 0.81° 0.86'' 



^ Combination of intersubband and interband spectroscopies and ab initio HSE06 + Go Wo calculations on 
GaN/AlN quantum wells with lattice-matched AIN from Ref. ^ J . 

ETBM + BZ average calculation, Ref. 32, no explicit specification on bulk parameters. 
Estimated after analysis of various experimental and theoretical sources, Ref. 

Difference of branch points from Ref. lit, see footnote of Tab. II. This should only be considered a rough estimate, 
as the branch points were originally obtained the opposite way around. 
° Difference of branch points from Ref. 17, see footnotes of Tab. II. 




Figure 3. (Color online) Band lineup for the zincblende 
modification of the Ill-nitrides. The small deviation from 
the intransitivity between the GaN/AlN, InN/GaN and the 
InN/AlN VBO originates from interface dipole effects (see 
text). 



lattice distortions is not of major importance when the 
electronic structure of the alloy is calculated. 

In fact, it has already been shown in Ref. that 
the CPA gives good agreement with supercell calcula- 
tions and further literature values for the F- valley bowing 
and the crossover concentration between a direct GaN- 
dominated and indirect AlN-dominated behavior. Be- 
sides a small variation in the band curvatures, which is 
not relevant for the here discussed properties, those cal- 
culations were performed for an almost identical param- 
eter set and the same VBO (albeit by chance, as the 
value of Ai;vB (GaN/AlN) sa 0.8 eV calculated in this 
paper coincides with the there employed recommenda- 
tion from Vurgaftman and Meyer). A closer look at the 
bandwidths and the position of the first band moments 
as depicted on the right of Fig. I shows that this system is 
very close to the weak scattering limit, which also under- 
lines the applicability of the CPA. Furthermore, recent 
experimental data for the band alignment is available. 



in particular for cubic nonpolar Alo.3Gao.7N/GaN het- 
erostructures, which were grown by means of molecular 
beam epitaxy on 3C-SiC (001). In Ref. 7, Wei et al. re- 
port on a two-dimensional electron gas at this interface 
and observe a large conduction-to- valence band offset ra- 
tio of 5:1 which enhances the electron accumulation. 

As the electronic structures of two materials A and B 
enter the calculation for alloys, an energy offset and thus 
VBO between A and B has also to be established as input 
parameter. In general, the choice of a VBO for CPA or 
supercell calculations for alloys is far from trivial. With 
varying concentration x in the A^Bi^x alloy, the out- 
line of the problem does in principle change from single- 
impurity-like behavior over dilute systems up to whole 
A-in-B-clusters, so that the "true" VBO should in prin- 
ciple be concentration- as well as lattice-site-dependent. 
Nevertheless, the usage of a fixed parameter is common 
for the sake of simplicity. In case of the AlGaN system 
under consideration, we will use the interface VBO of 
0.8 eV, as it has been shown to be directionally indepen- 
dent in the last section and gives the correct limit in the 
pure and in the phase separation case. 

The resulting CPA quasiparticle band structure for 
Alo.3Gao.7N is depicted in Fig. 4. Note that the 
linewidths were artificially enlarged for enhanced con- 
trast in this figure. For vanishing imaginary part in the 
complex energy z — E + iS, an unambiguous CB edge 
and VB edge and thus band gap can be identified. 

In order to obtain the position of the BP over the 
whole concentration range, a small imaginary part of 
6 = 10^* eV and #k w 10^ values in the irreducible 
wedge have been used in the CPA iterations, together 
with an energy resolution of Az = 5 x 10~^ eV. In the 
following, we restrict the discussion to the results of the 
BZ average approach as defined in Eq. (12) for two rea- 
sons: 

(1) The almost identical energy difference for pure 
GaN and AIN between Ebp{BZ average), i;BP (001) and 
-Ebp(IIO) of 0.1 and 0.25 eV, respectively, directly car- 
ries over to the results for AlGaN when the quasiparti- 
cle band structure is used according to Eq. (15). This 
holds in good approximation over the whole concentra- 
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Figure 4. (Color online) CPA quasiparticle band structure for 
Alo.3Gao.7N. 



tion range. Due to the very slow convergency, this has 
only been spot-checked with a lower energy resolution. 

(2) In real materials, disorder effects will prevent any 
translational invariance, which is an artificial symmetry 
of the CPA Hamiltonian. It is therefore at least question- 
able whether effects that strongly originate from symme- 
try properties of the quasiparticle band structure carry 
over to the behavior at real disordered interfaces. 

In addition, virtual crystal calculations have also been 
performed for comparison, as the VGA is also used in the 
literature to calculate branch-point energies of alloys, 
supposedly due to its simplicity. Here, the matrix ele- 
ments of the alloy with concentration x are obtained by 
a linear interpolation between the pure values. As the 
VGA is contained in the GPA as a limit case, the GPA 
results serve as a benchmark for the reliability of the 
VGA. 

The resulting band edges and BPs of AlGaN can be 
found in Fig. 5. Here, the same band subset as in the 
pure case was used for all concentrations, but for mate- 
rial combinations close to the split-band limit, an indi- 
vidual and concentration-dependent choice might be nec- 
essary. It is immediately evident that the position of the 
GPA as well as the VGA BP is almost constant on the 
common energy scale of the Aa;Bi_j, alloy, so that the 
relevant value for band alignment, the energetic distance 
between the BP and the VB edge, is only determined 
by the latter. This behavior originates in the use of a 
constant A-to-B VBO that has itself been determined by 
the alignment of BPs and will not hold when either a 
concentration-dependent VBO and/or a VBO with large 
dipole contribution is used in the calculations. 

A comparison of the GPA to the VGA results reveals 
that the latter gives slightly too low VB edges. In Fig. 6, 
the relative BP position E-Qpix) — E^j-Qix) is depicted as 
a function of the concentration x. While the VGA er- 
roneously yields an almost linear behavior, the relative 
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Figure 5. (Color online). CB edge, VB edge and BP energy 
of AlGaN, calculated in the CPA and the VGA from the TB 
band structures, respectively. The VGA results for the BP 
and the CB edge are overlaid by the CPA data. 
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Figure 6. (Color online) Relative BP position £'Bp(a::) — 
Ev-b{x) of AlGaN calculated in the CPA and the VGA. 



GPA BP position shows a pronounced bowing, where a 
satisfactory fit by a quadratic or cubic function is not 
possible for the whole concentration range. Table IV 
comprehensively lists the band edges and the relative 
BP energy of zb-AlGaN. These values can be used to 
estimate the band alignment for all concentrations with 
respect to any material whose branch points are known. 
Like in the pure case, the uncertainty will at least be 0.2 
eV, not including additional deviations due to the use of 
the GPA. 

We will finally compare our results to the work of Wei 
et al. for cubic Alo.3Gao.7N/GaN heterostructures. The 
difference in the electrostatic potential at the interface, 
which is determined by means of electron holography, 
gives a GB offset of 0.65 eV in their experiment. To- 
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Table IV. Concentration-dependent VB and CB edge -Evb/cb, 
and relative BP position Ebp — Evb for zb-AlGaN as obtained 
in the CPA. All values in eV. 



X \h\) 






-Ebp — -BvB 


1 





5.44^ 


3.28 


0.9 


0.20 


5.40 


3.09 


0.8 


0.30 


5.37 


2.99 


0.7 


0.39 


5.29 


2.91 


0.6 


0.46 


5.12 


2.83 


0.5 


0.53 


4.95 


2.77 


0.4 


0.59 


4.78 


2.71 


0.3 


0.65 


4.61 


2.65 


0.2 


0.70 


4.44 


2.61 


0.1 


0.75 


4.28 


2.55 





0.80 


4.11+ 


2.50 



Input values from parametrization, see Tab. I. 



and VBOs for alloyed systems and applied this method 
to the zincblende modification of AlGaN over the whole 
concentration range. While a virtual crystal treatment 
of the system gives erroneous results, the CPA agrees 
well with experimental data on Alo.3Gao.7N/GaN het- 
erostructures. Our results can be used to determine the 
band alignment in isovalent heterostructures involving 
pure cubic Ill-nitrides or AlGaN alloys for arbitrary con- 
centrations. 
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gether with a measured band gap difference of 0.78 eV 
across the junction, this yields a VB offset of 0.13 eV. 
The value is in very good agreement with our theoreti- 
cal result of 2.65 eV - 2.50 eV = 0.15 eV, However, this 
compliance definitely exceeds the predictive power of the 
calculations. When combined, our BP and band gap re- 
sults for this material combination give a conduction-to- 
valence-band ratio of 0.5/0.15 « 3.3, which is below the 
experimental value of 5. The deviation can be traced 
back to their spectroscopically determined band gap dif- 
ference between Alo.3Gao.7N and GaN, which is larger 
than our calculated band gap difference of 0.65 eV. Also, 
our calculations do not include any effects from lattice 
strain, which, albeit small, will also influence the results. 
Notably, the AlGaN layer thickness in their experiment is 
only 30 nm, so that finite size effects will also likely play a 
role. Furthermore, the authors notice composition fluctu- 
ations and oxygen contamination in their samples. If we 
account for these effects by using the same modified band 
gap difference of 0.78 eV in our calculations, we arrive at 
a ratio of 4.2, which is surprisingly good, especially when 
the various sources of uncertainty are considered on both 
the theoretical and the experimental side. 



IV. SUMMARY 

In this paper, we showed that the branch-point (BP) 
energies and resulting band alignment for the cubic III- 
nitrides InN, GaN and AIN can be calculated from tight- 
binding band structures with consistent and up-to-date 
input parameters. The directional independence of the 
valence band offset (VBO) at the unstrained GaN/AlN 
interface has been traced back to a constant shift of 
the corresponding BPs by the use of Tersoff's Green's 
function method. We then showed how the coherent 
potential approximation (CPA) can be used in combi- 
nation with the tight-binding model to obtain the BPs 
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